On the Bures Volume of Separable Quantum States * 

Deping Ye^ 

Abstract 

We obtain two sided estimates for the Bures volume of an arbitrary subset of 
the set of N x N density matrices, in terms of the Hilbert-Schmidt volume of that 
subset. For general subsets, our results are essentially optimal (for large N). As 
applications, we derive in particular nontrivial lower and upper bounds for the 
Bures volume of sets of separable states and for sets of states with positive partial 
transpose. 

PACS numbers: 02.40.Ft, 03.65.Db, 03.65.Ud, 03.67.Mn 

1 Introduction 

Quantum entanglement was discovered in 1930's [10, 28] and is now at the heart of 
quantum computation and quantum information. The key ingredients in quantum 
algorithms such as Shor's algorithm for integer factorization [29] or Deutsch-Jozsa 
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algorithm (see e.g. [24]), are entangled quantum states, i.e., those states which can not 
be represented as a mixture of tensor products of states on subsystems. Following [43], 
states that can be so represented are called separable states. Since determining whether 
a state is entangled or separable is in general a difficult problem [11], sufficient and/or 
necessary conditions for separability are very important in quantum computation and 
quantum information theory, and have been studied extensively in the literature (see e.g. 
[14, 15, 16, 17, 18, 19, 25]). One well-known tool is the Peres' positive partial transpose 
(PPT) criterion [25], that is, if a state on H = C Dl ® C° 2 ■ ■ ■ <g) C D " is separable then its 
partial transpose must be positive. Equivalently, if a state on 7i does not have positive 
partial transpose, it must be entangled. This criterion works perfectly, namely, the set 
of separable states S = S{7i) equals to the set of states with positive partial transpose 
VVT = VVT(H) for H = C 2 <g> C 2 (two-qubits), H = C 2 ® C 3 (qubit-qutrit), and 
H = C 3 <g> C 2 (qutrit-qubit) [14, 37, 44]. However, entangled states with positive partial 
transpose appear in the composite Hilbert space H = C 2 ®C 4 and TC = C 3 ®C 3 [15] (and 
of course in all "larger" composite spaces; see also [4] discussing the three-qubit case). 
One striking result is that, by some measures, the positive partial transpose criterion 
becomes less and less precise as N = Yl™ =1 V)% grows to infinity [2, 38]. This is inferred by 
comparing the Hilbert-Schmidt volumes of S and VVT , the estimates which rely on the 
special geometric properties of the Hilbert-Schmidt metric and were obtained by using 
tools of classical convexity, high dimensional probability, and geometry of Banach spaces. 
The same method can also be employed to derive tight estimates for the Hilbert-Schmidt 
volume of T> — ViTi) (the set of all states on 7i). However, a closed expression for the 
exact value of this volume is known; it was found in [45] via the random matrix theory 
and calculating some nontrivial multivariate integrals. 

2 



Compared with the Hilbert-Schmidt metric, the Bures metric on D [5, 42] is, 
in some measures, more natural and has attracted considerable attention (see e.g. 
[7, 8, 9, 20, 21, 39, 40, 41]). The Bures metric is Riemannian but not flat. It is monotone 
[26], i.e., it does not increase under the action of any completely positive, trace preserving 
maps. It induces the Bures measure [3, 12, 36], which has singularities on the boundary 
of T>. The Bures volume of T> has been calculated exactly in [36] and happens to be equal 
to the volume of an (N 2 — l)-dimensional hemisphere of radius ~ [3, 36]. (This mysterious 
fact does not seem to have a satisfactory explanation.) On the other hand, the precise 
Bures (or Hilbert-Schmidt) volumes of S and WT are rather difficult to calculate since 
the geometry of these sets is not very well understood and the relevant integrals seem 
quite intractable. These quantities can be used to measure the priori Bures probabilities 
of separability and of positive partial transpose within the set of all quantum states. 
(Here priori means that the state is selected randomly according to the Bures measure 
and no further information about it is available.) For small N, e.g., iV = 2 x 2 and 
N = 2 x 3, the Bures volume of S (hence of WT) has been extensively studied by 
numerical methods in [30, 31, 32, 33, 34, 35]. For large N, the asymptotic behavior of 
the Hilbert-Schmidt volume of S and WT was successfully studied in [2, 38]. Based 
on that work, we shall derive in this paper qualitatively similar "large N" results for 
the Bures volume. In summary, our results state that the relative size of S within T> is 
extremely small for large N (see Corollaries 1 and 2 for detail). On the other hand, the 
corresponding relative size for WT within T> is, in the Bures volume radius sense (see 
section 2 for a precise definition), bounded from below by a universal (independent of N) 
positive constant (see Corollary 3). The conclusion is that when N is large, the priori 
Bures probability of finding a separable state within WT is exceedingly small. In other 
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words, we have shown that, as a tool to detect separability, the positive partial transpose 
criterion for large N is not precise in the priori Bures probability sense. Its effectiveness 
to detect entanglement is less clear (see the comments following Corollary 3). 

This paper is organized as follows. In section 2, we review some necessary 
mathematical background, particularly the background for the Hilbert-Schmidt volume 
and the Bures volume. Precise statements of our main results can be found in section 
3. Section 4 explains why our estimates are essentially optimal for general subsets of 
quantum states. Section 5 contains conclusions, comments and final remarks. 

2 Notation and Mathematical Preliminaries 

2.1 Mathematical framework 

We now introduce the mathematical framework and some notation. Let 7i be the 
(complex) Hilbert space C Dl ®C D2 • ■ -®C D ™ with (complex) dimension N = DiD 2 ■ ■ ■ D n . 
Here we always assume n > 2 and Di > 2 for all i = 1, 2, • • • , n. Recall that Di = 2 for all 
i corresponds to n-qubits, and Di = 3 for all i corresponds to n-qutrits. n = 2 corresponds 
to bipartite quantum systems and n > 2 corresponds to multipartite quantum systems. 
Denote by B(?i) the space of linear maps on 7i. Define the Hilbert-Schmidt inner product 
on space B(TC) as (A, B)hs = tr(^5). The subspace oiB(Tl) consisting of all self-adjoint 
operators is B sa (Ti). It inherits a (real) Euclidean structure from the scalar product 
(■,-)hs on B(H). (This is because if A, B e B sa (H), then (A,B)hs must be a real 
number.) V denotes the set of all states on 7i (more precisely, states on B(TC)), i.e., 
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positive (semi) definite trace one operator in £> sa (7i): 



V = V(H) := {p £ B sa (H),p> 0,trp = 1}. 



A state in D is said to be separable if it is a convex combination of tensor products of n 
states (otherwise, it is called entangled). Denote the set of separable states by S, then 



Both V and S are convex subsets of B sa (H) of (real) dimension d = N 2 — 1. 

Indent: The notation S(7i) is in principle ambiguous: separability of a state on 13(H) 
is not an intrinsic property of the Hilbert space TL nor of the algebra B(TC); it depends 
on the particular decomposition of Ti, as a tensor product of (smaller) Hilbert spaces. 
However, this will not be an issue here since our study focuses on fixed decompositions. 

2.2 Hilbert-Schmidt and Bures Measures on V 

Any quantum state on Ti can be represented as a density matrix, i.e., the N x N positive 
(semi) definite matrix whose diagonal elements sum up to 1. Therefore, any quantum 
state p eD has eigenvalue decomposition p = UAW for some unitary matrix U G U(N) 
and some diagonal matrix A = diag(Ai, • • • , \n) with (Ai, • • • , \n) £ A. Hereafter, Idjv 
is the N x N identity matrix and U £ U(N) means that U is an N x N matrix with 
UIP = U^U = Id at- We denote by A the regular simplex in R N , i.e., 



The Weyl chamber of A defined by the constraint Ai > ■ ■ ■ > A^ is denoted by Ai. 
Clearly, for any p = UAW as above and for any diagonal matrix B £ U(N), we 



S = S(H) := conv{p! <g> • • • ® p n , pi £ V(C Di )}. 



A 
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have UAU^ = UBAB'W . Thus, to have unique parametrization of generic states 
p = UAW G T>, we have to restrict (Ai, • • • , Ajy), for instance, to Ai and select one 
specific point in the coset space F N = U (N) / [U (1)] N (the flag manifold). 

We will be interested in various measures on P. A natural restriction is to require 
invariance with respect to unitary rotations. For most problems, the interesting class 
of measures are those that are invariant under conjugation by a unitary matrix. Such 
measures can normally be represented as the product of some measure on Ai and the 
invariant measure on T N (see [3, 13] for more on this and for the background on the 
discussion that follows). The unique (up to a multiplicative constant) invariant measure 
7 on T N is induced by the Haar measure on the unitary group U(N) and has the form 

d 7 = 2Re{U- 1 dU) ij Im(U- 1 dU) ij , 

X<i<j<N 

where U G U(N) and dU is the variation of U such that U + dU G U(N). The total 7 
measure of T N is known to be (see [45]) 



(27T) 



N(N-l)/2 



Zn= K ' , where E(N) = l% =1 T(j). (1) 

Here T(x) = / °° t x ~ 1 e~ t dt is the Gamma function. 

The Hilbert-Schmidt measure Vhs(-) on V, induced by the Hilbert-Schmidt metric, 
may be expressed as [45] 

N-X 

dV HS = VN J] (Aj - Aj) 2 Y[ d\ % d 1 , (2) 

X<i<j<N i=l 

where (Ai, • ■ ■ ,Xn) G Ai. (This is just a different name for the canonical <i-dimensional 
Lebesgue measure on D.) Therefore, to obtain the Hilbert-Schmidt volume of T>, one has 
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to calculate the following integral [45]: 



» JV-l 

V HS (V) = / v/iV (A, -A,) 2 ] 

J A! x F» l<i<j<N i=l 



We define vrad#s(/C), the Hilbert-Schmidt volume radius of JC C T>, to be the radius of 
(i-dimensional Euclidean ball which has the same volume as the Hilbert-Schmidt volume 
of /C. In other words, 

vrad# s (/C) = 



where ad = T (J^/2) * s ^ ne vomme of d = N 2 — 1 dimensional Euclidean ball. For later 
convenience, we also denote VR H s(JC,C) as VR H s(JC,C) = ( vgs(ic) ) = vradggfS ' lt 
amounts to comparing the Hilbert-Schmidt volume radii of JC and C. 

It is known that vrad#s(2}) ~ e~*d~^ [38] by Stirling approximation 

rW = #mi + ofiY). (4) 



z \e/ \ \ z 



Stirling approximation (4) also implies that (ad) d ~ \/2en d 2 and therefore 

~ (47r 2 e)3d-f. (5) 

Here a(n) ~ b{n) means lim n ^ 00 a{n)/b{n) = 1. 

An arguably more important measure in the present context is the Bures measure (or 
Bures volume) Vb(-), which can be written as [36] 

dv B = 2 n 4-^11^^ (6) 

VAi • • • Ajv a * + a j i=1 
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where (Ai, • • • , \n) G Ai. The Bures measure is induced by the Bures distance ds( 



which may be defined via ds{Qi, Q2) — y 2 — 2tr^/ ' y/~QiQ2y/Qi f° r an y states Qi, Q2 € D. 

The Bures measure has singularities (with respect to the Hilbert-Schmidt measure) 
on the boundary of T>. (The boundary corresponds to at least one of the Aj's being 0, 
and if two or more of them are 0, then some denominators in (6) vanish.) Thanks to the 
work of Sommers and Zyczkowski [36] , we know the precise value of the Bures volume of 
T>, that is 

v B (v) = / j=== n %^n^ 



,1-AT 2 ^ 



(7) 



r(iv 2 /2)' 

As mentioned earlier, this value happens to be the rf-dimensional volume of the d- 
dimensional hemisphere with radius \. We define vrad#(/C), the Bures volume radius 
of /C C V, to be 

vradB(;c) = (M2)\ 

While comparing the Bures volume of /C with the Hilbert-Schmidt volume of the 
Euclidean ball does not have immediate geometric meaning, we find this way of describing 
the size of K, in the Bures volume sense convenient in our calculations. 

By formulas (4) and (7), one has vrad^'C) ~ | and hence 

(V B (V))^^ d -h. (8) 

For later convenience, we also define the (relative) Bures volume radii ratio of /C to £ as 

VRb(/C, C) = ( yg(£) ) — vrachoc) ' can ^ e use< ^ as a measure °f the relative size 
of /C to £ in the Bures volume sense, and clearly does have geometric meaning. 
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We refer the reader to the references [3, 12, 13, 36, 45] for more detailed background 
and for motivation. In the following sections, we are interested in the (asymptotical) 
behavior of VRb(/C, V) in terms of its relative VRhs^/C, V). 



3 Main Results 



In this section, /C will be an arbitrary (Borel) subset of T>. We will estimate the Bures 
volume of /C, in terms of the Hilbert-Schmidt volume of /C, both from below and from 
above. The following lemma is our main tool to study the asymptotical behavior of 
VRb(/C,P). We point out that these estimates are independent of the possible tensor 
product structure of 7i. 



■^V B {K)<{Y^Yl{p)-^. 



Lemma 1 For any subset K, inT> and any p > 1, one has 

JV-JV 2 AT 2 -1 , J 

2S- Ns-V H s(1C) < 2" 
where I(p) is defined as 

W- m r (tap) 111 r(|5f) ) EW ■ U 

Remark. I(p) can be defined for all p ^ [|, 1] (irrespective of N; since the Gamma 
function has poles at nonpositive integers, there are singularities in [|, 1] whose exact 
locations depend on N.) In particular, 1(0) = ^^=2^. The quantity E(N) was defined 
in (1). 

Proof. First of all, we estimate Vb(JC) from below. To that end, define h : A — > R as 

N 

t=l l<i<j<N 
9 



Lagrange multiplier method implies that (1/iV, ••• ,1/N) is the only critical point of 
h(Xi, • • • , Xn) in the interior of simplex A. Clearly h(Xi, • • • , Xn) is always on the 
boundary of the simplex A, which consists of sequences for which one or more of the 
Aj's equal to 0, and strictly positive in the interior of the simplex A. By compactness, 
h(X\, • • • , Xn) must have a maximum inside, and the critical point (1/iV, • • • , 1/N) must 
be the (only) maximizer of h(X\, • • • , Xn) on A. Therefore, 

1 1 -I— r 1 N-N 2 N 2 

II rzr^ 2 " 5- ^ ( 10 ) 



y/h(Xi,--- ,X N ) VXi---X N ^ i<:j ^ N \ + \ 
By formula (6), the Bures volume of /C equals to L. cLVb, i.e., 

Jk v A i--- A ^i<^i<jv Ai + A ^ tl 

Considering inequality (10) and formula (2), one gets 

2 r N ~ x 
V B {K) > 2 X - N2 N- J] (Xi - Xj) 2 dXi dj 

•* K Ki<j<N i=l 



2^ N2 N^V HS {K). 



Next, we will derive the upper bound, which is more involved (and more important 
for our results). The subset <9A, the boundary of A, consists of sequences for which some 
Xi = and has zero N — 1 dimensional measure. Thus, without loss of generality, we 



can assume Aj > for all % — 1, • • • ,N and, in particular, 



< 1 for all i ^ j. This 
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implies 

N 2 +N _ 2 r l „ (Aj - Xj) 



2 N-l 



Ik VAi • • • X N Ai + Aj f£ 

AT-1 



< 



/AC \/Ai • • ■ Aat 



|Aj - Aj| ] [ dA^c^y 

l<i<j<N i=l 



AT-l 



'AC 



/ 5r JJ dA^ 7 (11) 

- i=i 

where, to reduce the clutter, we denoted 



ff(Ai, • • • , Ajy) = — == (Ai-Ajl 1 p. 

V Ai • • • Ajy 

/(Ai, • • • , Ajy) = J] l^- A il"- 

l<i<j<N 



For any p > - (so that 2p > 1), we employ the Holder inequality to (11) and get 

2 iL ^ 1 V B (IC)< [JJ 2P TI *>«*l) P [J k 9^U dXid^f) " . (12) 
Substituting / into the first integral of (12) and by (2), one has 

(// 2 *ff^)*=(7 n (A,-A/n^f=(%Bf. 

\JlC i=1 J ^ JlC l<i<j<N i=l ' ^ / 

Substituting g into the second integral of (12) leads to 

/JV-l „ . JV-l 

i=l l<i<j<N i=l i=l 

/JV p _ 1 JV-l 

l<i<j<N i=l i=l 

the inequality following just from K <Z T>. By (1) and the Fubini's theorem, the last 
integral in (14) equals to 

^ Al l<j<j<iV i=l i=l 



£(iV) 
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Under the condition ^—^ > (i.e., p > 1 or p < 1/2), one has (see e.g. [22, 45]) 



I n i^-i-n# i " I n^=- 7 ^ T (n 

- /A l<i<j<7V <=1 i=l 1 ^ 2p _! j S=l 



r-r 1 l^ 1 ^ 2p-l y 1 ^ 2p-l 



r ( 3p ~ 2 

1 2p-l 



Taking into account that A consists of AM Weyl chambers, we conclude that the expression 
in (15) is then equal to I(jp). In other words, we have shown that 

„ TV _ 1 2V-1 

i(p)= n i^-^-i^n^^n 

l<i<j<N i=l i=l 

Combining this with (12), (13), and (14), we conclude that if p > 1, then 

2^^(K) < (^^)%)^, 
which is the upper estimate from Lemma 1. 

Theorem 1 There is a universal computable constant c\ > 0, such that for any Hilbert 
space Tt and any subset JC CD, 

c x VR HS (IC,V)<VR B (IC,V). 

Proof. Recall d = N 2 — 1. From the lower bound of Lemma 1, one has 

V B {K,f* > l - V HS {K)* >\d^ Vhs{K)*. 

Dividing both sides by V B (V)^, one obtains 

V B (JQ\* > l d i (V HS {K,)\* (V HS {V)\^ 



V B {V)J 2 \V HS (V)J \V B (V) 

1 i / vrad H5 (P) \ / Vg g (/C) \ 3 

2 4 V vrad B (X>) J \Vhs(V)J 



(16) 
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Formulas (5) and (8) imply that d* vrad#s(r>) ~ 2e 4vrad.e(X>), i.e., 



lim - c? 4 

d— >oo 2 



vrad gg (T') 
vrads(P) 



) 



= e 4. 



Therefore, there is a (computable) universal constant c\ > 0, such that, for any N, 



compute ci precisely. For instance, ci(4) ~ 0.7572. Actually, ci(6) ~ 0.7686, ci(8) ~ 
0.7728 and Ci(10) ~ 0.7748 which are very close to the e^ 1 / 4 ~ 0.7788. It appears that 
the sequence C\(N) is increasing and so c\ = ci(4) should work for all N, however, we do 
not have a rigorous proof. 

Theorem 2 There is a universal computable constant C\ > such that, for any Hilbert 
space H and any JC C V, 



where a = VR#s(/C, P). 

Proof. Recall iV = DiD 2 ■ ■ ■ D n and d = N 2 — 1. For any p > 1, Lemma 1 implies that 



1 i vrad g g(£>) 

2 vrad B ('D) 



> ex. 



Together with (16), this shows that for any /C C V and for any iV, 



VR B (}C,V)>c 1 VR HS ()C,V). 





d — 2 d 1 2p— 1 

22Vb(/C) < 2 — a — Vb(/C) < {V HS {V))^ I(p) — . 
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Let (3 = Replacing Vhs(^) and I(p) by formula (3) and formula (9), one has 



2*V B {K) < (2n 



[AT! 



[r(A^ 2 )] 2? [r(/3iv 2 )r(i + /3) 



Clearly [^(A^)]^ 1 < 1 and N^(Nl)^' 1 < 1 if p > 1. Hence, 

[nL (r^ra + A?))] 1 ^ 



Vb(/C) < a 2 ? 7T2 



[r(iv 2 )]5? [r(/?jv 2 )r(i + /?)"] 



(17) 



2;> 



Since xT(x) = F(x + 1), it is easy to see that for all x G (0, 1) the upper estimate 
of T(x) is - and (somewhat less easy that) the lower estimate of T(x) is 4-, where 
•d « 1.12917 [6]. That is 



— <r(x)<~, or -<r(l + x)<l, for all x G (0, 1) 

VX XV 
N 2 ln(e/q)-l 



'11 



Pick p = p(N, a) := ^2 i n ( e / Q j_ 2 as a function of N and a, so that (3 = ^ = N 2 ln{e/a y 
Equivalently, (5N 2 = ln(e 1 /a) . Since a < 1 for all fC CV, then /3iV 2 < 1 and hence /3j < 1 
for all j = 1, 2, • • • iV 2 . Taking inequality (18) into account, one has 



r(i + A?) 



r(i + /?) 

Consequently, again by inequality (18), and N > 4 



for aU j = 1,2, ■■■N. 



(19) 



T(pN 2 ) 



2p 



< 



■& N 



{N-l)\ 



2p 



p(i-m-& < /3 (i-jv)(i-^) > (20) 



Combining inequality (17) with inequalities (19) and (20), one has 



V B {K) < 7T2 iT (3 



N fl (l-JV)(l-i) 



[r(7v 2 )] 



Equivalently, taking d-th root from both sides, 



(Vb(£))3 < ^ [r(iv 2 )]^, 
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Note N > 4, and hence fin-i < $3 = 1.0413. Now dividing VB{T>)d from both sides of 
the above inequality, one gets 

VR B {JC,V) < 2^3 a& [iV 2 ln(e/a)]^r (1 ~^ ) [r (N 2 /2)}^ [T(N 2 )}^ . (21) 

It is easy to verify that 

rAT o, , , N1 1 /lniV\ /lnln(e/a)\ r /lnln(e/a)\ 

[iV 2 ln(e/a)] < exp l^—) exp ( ^ < V2 exp f ^ } j , 

Therefore, 

[N'He/a)}^ 1 -^ <V2e^ exp (^^) (22) 

Also, we can verify that 

jl 1 / ln(l/a) \ 1 / 1 \ jl r- , s 

a2p=a ~ 2 exp U^(i + HIM) -i] )- a ~ 2 exp \w) ~ 632 ^ (23) 

Since r(iV 2 ) < (iV 2 )! < exp(iV 2 In iV 2 ), one has 

ln(r(iV 2 )) 



1 1 



T(N ) 2d 2 P d = exp 



2(N 2 - 1) (iV 2 ln(e/a) - 1) 
21niV \ / 2 , 

£ exp <iVM^J £exp U> £e '' (24 > 

Stirling approximation formula (4) implies that 

lim -i '—L- = —. 25 

n^oo [r(N 2 )}^ V2 

Together with inequalities (21), (22), (23), and (24), there exists a universal (independent 



of N, a) constant C\ > 0, such that, VRb(/C, X>) < C\y/a exp ( 



In ln(e/a) 
2JV 
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Remark. A slightly more precise calculation shows that 



VR B (£,£>) < V2a exp 



In ln(e/a) 
2A 



1 + 



In A 



The calculation yields explicit (not necessarily optimal) values of C\ in the theorem. 
For small dimensions, our proof yields Ci(4) ~ 2.5164 if A = 4, Ci(6) ~ 2.2137, and 
Ci(8) ~ 2.0478. As the dimension A becomes large, the value of C\ given by the argument 
tends to y/2 ~ 1.4142. On the other hand, the Legendre duplication formula (see [1]) 
says that 

r(z) r (z + 1/2) = 2 x - 2z v^F T(2z). 
By taking z = N 2 /2, one can rewrite the expression in (25) as 



r(A 2 /2) r(N 2 /2)\^^) 



ttT(N 2 /2) 



T(A 2 ) J ^2 Vr(A 2 /2 + l/2), 

Gamma function is log-convex [23] , and hence 

T(A 2 /2 + 1/2) 2 



(26) 



T(A 2 /2) 2 < T(A 2 /2 - l/2)r(A 2 /2 + 1/2) 



A 2 /2 - 1/2 



Equivalently 



T(A 2 /2 + 1/2) 



> 



A 2 



r(A 2 /2) - V 2 

1 « 2.7. Together with formula (26), this shows 



which is greater than y/n iff A > \/2tt 
that the asymptotic relation (25) is in fact an upper bound for all A > 3. It follows that 
d w 2.5164 works for all A > 4, C x » 2.2137 works for all A > 6, etc. 

Remark. In most cases of interest a is such that the factor exp y ^^ j is bounded 
by a universal numerical constant. For instance, if ln(l/a) < a% e a2N for some constants 
ai > 0, a 2 > 0, then 

'In A 



VR B {)C,V) < V2^ y/VR HS (}C,V) 



1+0 



N 
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While our argument doesn't give similar estimates for general a, other ways of writing 
the estimates in more transparent ways are possible. For example, for any fixed p > 1 
there is a constant C p > depending on p (but independent of N and a) , such that 

VR B (/C,P) < C p (YR HS ()C,V))^. 
In the cases of S and WT, - is bounded from above by N k for some (fixed) integer k 



[2, 38]. Therefore, VRb(«S, T>) < C\ a/VRh-s(iS, V) where C\ > is a universal constant 



independent of N. Similarly, VR B (VVT,V) < C x ^/VR H s(VVT,V) where d > is a 
universal constant independent of N. 

Remark. We point out that there is a lot of flexibility in the choice of j3 — N $ Yn(e/a) 
(hence the choice of p(N,a)). For example, one can choose (3 = eN ^?r^ , and proves 
Theorem 2 with different (larger) constants. However, formula (23) does suggest that the 
factor ln(e/a) in j3 is essentially optimal in general. 

As applications of Theorems 1 and 2, and the estimates for VRhs{<S, V) implicit in 
[2], one immediately has the following corollaries. 

Corollary 1 (Large number of small subsystems) For system 7i = (C D )® n , there 
exist universal computable constants C2, C2 > 0, such that for all D,n > 2, 



<VR B (S,V)<C 2 -l {Dnlnn)1/2 



where a D = \ log D (l + -5) - jm 1o Sd{ d + !)• 

Corollary 2 (Small number of large subsystems) For system Ti = (C D )® n , there 
exist universal computable constants C3, C3 > ; such that for all D,n > 2, 



c? TTT ^ ^ „ \ (nlnn) 1 / 2 

< VR B (S,V) <C 3 ■ 



N l/2-l/(2n) ~ «V ? ; _ ^ y jyi/2-l/(2n) ' 
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Remark. Recall that if 7i = (C D )® n , then the dimension of H is iV = D n and so 
the expressions in the the numerators of the estimates in the Corollaries above are of 
smaller order than the denominators. Hence, for any fixed small D, Corollary 1 shows 
that VRb(iS,P) goes to exponentially asn^ oo. On the other hand, for H = (C 13 )®™ 
and for fixed small n, Corollary 2 shows that "the order of decay" of VRb(S,T>) is 
between _D^~f and as D -> oo. In both cases, the priori Bures probability 

of separability is extremely small for large (and even for moderate) N. It is possible to 
provide (not necessarily optimal) estimates on the constants appearing in both corollaries. 
For instance, in Corollary 1, one can take c 2 (4) = 0.2272, C 2 (4) = y/4A d{4) = 5.2785 
if N = 4, c 2 (6) = 0.2306, C 2 (6) = V^A Ci(6) = 4.6436, and c 2 (8) = 0.2318, C 2 (8) = 
\/AA Ci(8) = 4.2955. If the dimension N is large (particularly for large n), the relevant 
asymptotic behaviors of c 2 and C 2 given by the proofs are: c 2 tends to ~ 0.6577, and 
C 2 tends to 2 3 / 4 e 1/8 « 1.9057. In Corollary 2, c 3 can be taken as c 3 = e~ 1/4 /VQ ~ 0.3179 
and C3 = C 2 . We refer the readers to [2] for the constants in terms of the Hilbert-Schmidt 
volume. 

In the rest of this section, we will discuss the Bures volume of WT . For a bipartite 
system 7i = C Dl <S> C D ' 2 , any state p on Ti can be expressed uniquely as 

Dl D 2 

where {e^}^ and {f a }a=i are the canonical bases of C^ 1 and C D2 respectively. Define 
the partial transpose T(p) with respect to the first subsystem as 

T (P) = ^2PjaM e i ® fa){ej ® //?!■ 
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We write VVT for the set of states p such that T{p) is also positive. (Note that VVT 
is basis-independent because eigenvalues do not depend on a basis [25].) The Peres 
criterion asserts: every separable state has a positive partial transpose [25]. That is, 
S C VVT C V. For qubit-qubit system C 2 ® C 2 and qubit-qutrit (or qutrit-qubit) 
system C 2 ® C 3 (or C 3 <E> C 2 ), the positive partial transpose criterion gives a sufficient 
and necessary condition for separability, i.e., S = VVT [14, 37, 44]. The following 
corollary, which is a direct consequence of Theorems 1, 2, and of Theorem 4 of [2], gives 
the estimation of the Bures volume of VVT . 

Corollary 3 (Bures Volume of VVT): There exists an absolute computable constant 
Co > 0, such that, for any bipartite system 7i = C D ® C D , c < VRb(VVT,V) < 1. 

Indent: An unsolved question is: does there exist a universal constant < Cq < 1 such 
that VRb(VVT,V) < Co < 1 ? Answering this question would help us understand 
the effectiveness of positive partial transpose criterion as a tool to detect quantum 
entanglement for all D > 3. (The answer to the analogous question about VRhs is 
not known, either.) 

An immediate consequence of Corollaries 2 and 3 is that, for Ti, = C D ® C D and 
large D, there exist universal constants c 4 ,C 4 > (independent of D), such that, 
c±D-2 < VR B (S,VVT) < C A D~4. The upper bound decreases to as D — > 00. In 
other word, the conditional priori Bures probability of separability given positive partial 
transpose condition is exceedingly small. Hence, for large N, the PPT criterion is not 
precise as a tool to detect separability. 
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4 Optimality of the bounds 

In this section, we will prove that, in general, the bounds in Theorems 1 and 2 are 
essentially optimal. The Bures volume has singularities close to the boundary of T>, so 
the optimal upper bound is intuitively attained by the subsets close to the boundary of 
V. On the other hand, for the subsets located near the maximal state p ma x — we 
can achieve the lower bound (this is really a simple consequence of the proof of Lemma 

!)• 

4.1 Optimality of the lower bound 

For < t < 1, let K t = W + (1 - t)p maxj i.e., 

K t = : X = diag H-^- + t\ x , • • • , ^ + t\ N ^j , (A l7 • • • A at) 6 A and U G W(iV) J 

Let Ztv be as in (1). 

We now estimate Vb(^Q) from above. By formula (6), 

N*+N-2 /" Zat t^ 2 " 1 -p-r (Aj-Aj) 2 TT , x 



As ^ + tAj > for all i, one obtains 



2^-V B ()Q) < Z N 2(N2 _ N)/2 (1 _ f)N2/2 J] (A 4 -A,) 2 n^ 



l<i<j<N i=l 



,7V2_1 N (N 2 -l)/2 

V HS (D), 



2(N 2 -N)/2 n _ f^N 2 /2 

where the equality follows the formula (2). By Lemma 1, one gets 

jJV 2 -l j.N 2 -l 
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Hence, VR B (JC t ,V) < ^ < At for all t < §. On the other hand, VR H s(JC t ,V) = t holds 
trivially because of the homogeneity of the Hilbert-Schmidt measure. We have proved 
that VR B (JC t ,V) < 4VR HS (JC t ,V) for all fC t such that VR H s(IC t ,V) < §. 

Theorem 1 guarantees that the lower bound of VR B (fC t ,T>) is at least (up to a 
multiplicative constant) VRf/s^/Q, £>). So the lower bound of "VTIb(/Q,:D) in Theorem 1 
can be obtained, and hence is optimal in general. 

4.2 Optimality of the upper bound 

For < t < 1, we consider K} as 

K} = {UXU ] : X = diag(l - t + t\ 1 ,t\ 2 , • • • , Ujv), (A 1; • • • X N ) e A 1 and C/ e W(iV)}. 

Recall Ai is the chamber of A with order Ai > ■ • ■ > A at. 

The Hilbert-Schmidt volume of JO can be calculated by the following integral 

N N-l 

VnsiJC 1 ) = Z N t^ 2 VN J] (A, - Xj) 2 HitX, - tX k + 1 - tf J] dA,. 

^ Al 2<i<j<N k=2 i=l 

As < tAi - tX k + 1 - t < 1, one has 




< ^2 ^(l"-" 2 V„ s (»«-i). (27) 

— 2 

By Stirling approximation (4), VR#5(/C*, £>) < C 4 t t'v+r holds for some universal 
constant C4 > 0. If t N + 1 is bounded from above, e.g., t > e C4 <~ 1_Ar ) for some constant 
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C4 > 0, then 

YR HS {K\V)<C 5 t (28 
holds for a new universal constant C5 > 0. 

Next, we estimate the Bures volume of /C* from below. By formula (6), 

jV 2 + AT-2 , + . 

(iv-D 2 2 N , . , x 2 AT-1 

11 x. j_ x. 11 Z7 11 rfA * 



Al v/^Ai + l- ^Aa-.-Ajv 2 <4<jv Ai + ^ fc= 2 tAl + tXk + 1 " * 

y Al v A2 • ■ ■ \n 2 <t4<N Ai + A ^' r=i 

where the inequality is because of < t\i + 1— t < 1, < tX\ + tXk + 1 — t < 1 and 
tXi — tXk + 1 — t > 1 — t. The last integral can be computed as in (27) and leads to 

N 2 +N-2 . t . 

2—^- V B {K}) 

Zm , , », r, (JV-1) 2 jv 2 _jv_2 s , , iV 2 -2]V-l 

>-^(l-t) 27V " 2 t^2^ 2 — Vk(^-i) / (l-Ai)^s— dAi 
Av-i 







^-1 (iv - l) 2 

Employing the Stirling approximation (4) one gets VRb(/C*, T>) > c$ \/t for some universal 
constant C5 > if, say, t < Together with (28), we have thus proved 



VR B (/C',P) >C5VVfW/C*,D) 
for some universal constant C5 > 0, if t < I and t > g^C-i-w^ Theorem 2 guarantees 



that VR B (/C*,X>) < Ci V VR H5(/C 4 ,P). Therefore, the upper bound in Theorem 2 can 
also be achieved, and is optimal in general. 
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5 Conclusion and Comments 



In summary, we proved that if K, is a Borel subset of V, then the priori Bures probability 
of /C can be estimated from above and from below in terms of the priori Hilbert-Schmidt 
probability of /C. Specifically, under some mild conditions on K, the relative Bures volume 
radius VRb(X,P) can be (approximately) bounded from below by the relative Hilbert- 
Schmidt volume radius VR#,s(/C, and from above by ^/VR#s(/C, V). We employ 
these results to estimate the Bures volume of S and VVT and the relevant priori Bures 
probabilities. We deduce that positive partial transpose criterion becomes less and less 
precise as the dimension of Ti becomes larger and larger, at least if the goal is to detect 
separability. We also give examples showing that, for general subsets, our bounds are 
essentially optimal. 

When N is small, for instance when N = 4 or 6, our estimates for Vb{S) are less precise 
than Slater's numerical results. However, our methods overcome the big disadvantage 
of the numerical approach, which works only for small N. Moreover, our results are 
independent of the structure of Ti. (In applications, of course, the information on the 
structure of Ti will be hidden in the calculation of VR/fs(/C, V).) Proceeding along similar 
lines one can obtain similar results for real Hilbert spaces, and then estimate the Bures 
volume of S or VVT on a real Hilbert space. 

As is well known, for sets in a Euclidean space (in particular, for sets of matrices 
endowed with the Hilbert-Schmidt metric) the volumetric information is roughly 
equivalent to the metric entropy information such as covering and packing numbers 
(see, e.g., [27]). However, for the Bures geometry the parallels are not so immediate. 
Consequently, further work is required to answer (even approximately) questions of the 



23 



type: Given e > 0, what is the maximal cardinality of a subset of T> (or S), every two 
elements of which are at least e apart in the Bures metric? 
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